Infrared Gluon and Ghost Propagator Exponents From Lattice QCD 



00 

o 
o 

(N 

Or 
<D 
00 

^1" 

(N 



> 

q 

in 
o 

o 



0. Oliveira, P. J. Silva 
Dep. Fisica, Universidade de Coimbra, 3004-516 Coimbra, Portugal 
(Dated: September 24, 2008) 

The compatibility of the pure power law infrared solution of QCD and lattice data for the gluon 
and ghost propagators in Landau gauge is discussed. For the gluon propagator, the lattice data is 
well described by a pure power law with an infrared exponent k ~ 0.53, in the Dyson-Schwinger 
notation, k is measured using a technique that suppresses finite volume effects. This value implies a 
vanishing zero momentum gluon propagator, in agreement with the Gribov-Zwanziger confinement 
scenario. For the ghost propagator, the lattice data seem not to follow a pure power law, at least 
for the range of momenta accessed in our simulation. 

PACS numbers: 12.38.-t, 11. 15. Ha, 12. 38. Aw, 14.70.Dj 



I. INTRODUCTION AND MOTIVATION 



The infrared properties of the Landau gauge gluon and 
jhost propagators in momentum space, respectively, 



G ab (q) = -S ab G(q 2 



8, v - ) D(q 2 ^ 



(1) 
(2) 



are connected with gluon confinement mechanisms, 
namely the Kugo-Ojima scenario (KO) [1] and the 
Gribov-Zwanziger horizon condition (GZ) The GZ 

mechanism requires D(0) = (which implies maximal 
violation of reflection positivity) and an enhanced ghost 
propagator, relative to the perturbative function. The 
KO confinement mechanism demands l/(q 2 G(q 2 )) = 
in the limit q — > 0. From the point of view of the KO 
and GZ confinement mechanisms, the requirements on 
-D(O) and G(0) are necessary conditions and its violation 
immcdiatly rules out these scenarios. 

In the recent years there has been a renewed interest 
in the computation of gluon and ghost propagators in 
the pure gauge theory, due to progress on solutions of 
the Dyson-Schwinger equations (DSE) and lattice simu- 
lations which explore further the infrared region. 

The DSE are an infinite tower of coupled nonlinear 
equations for the QCD Green's functions. The compu- 
tation of the gluon and ghost propagators within DSE 
requires denning a truncation scheme and parametrizing 
some of the QCD vertices. For the Landau gauge, recent 
solutions can be classified in two categories: solutions 
which are not compatible with the KO or GZ confining 
mechanism [J, |5j, |g, ljj and solutions which do not rule out 
the two confining mechanisms 

BSE, EH- For the first 
class of solutions, D(0) is finite and does not vanish. For 
the later class of solutions, in [l^] an analytical solution 
was found for the deep infrared region [l3| . The solution 
assumes infrared ghost dominance and connects the two 



propagators via a single exponent, n, 



Z(q') = q'D{<f) = CJ \ 



F(q 2 ) = q 2 G(q 2 ) = J \ 



(3) 
(4) 



a is a constant with dimension of mass. Moreover, DSE 
equations predict k — 0.595, which, for the zero momen- 
tum, implies a null (infinite) gluon (ghost) propagator. 
Furthermore, renormalization group analysis [141 |l5j re- 
strict the possible values for k to 0.52 < n < 0.595. 
This result suggests a null (infinite) zero momentum 
gluon (ghost) propagator. In [TtJ it was argued that 
the solution ©-(U]) is the unique power law infrared so- 
lution compatible with DSE and functional renormal- 
ization group equations. In (l6j a similar analysis of 
the DSE within time- independent stochastic quantisation 
predicted the same behaviour and k = 0.52145. 

The computer simulations of 4D pure SU(3) and SU(2) 
gauge theories on a lattice HJJlljM III III III III lM 
111, S3, H, US H [Ml, [H, m, \M, llj show a finite non- 
vanishing zero momentum gluon propagator. Although 
D(0) is a decreasing function of the lattice volume, from 
naive extrapolations to the infinite volume [29T I32I [33j it 
is not clear if the lattice result becomes compatible with 
the GZ confinement mechanism, i.e. if -D(O) approaches 
zero as V — > 00. Moreover, recent 4D lattice QCD cal- 
culations of the gluon propagator on large volumes, (27 
fm) 4 for SU(2) [H and (13 fm) 4 for SU(3) [28], show 
a propagator which goes to a constant in the infrared 
limit, although not excluding a gluon propagator tend- 
ing towards zero. However, to work out such huge lattice 
volumes, the simulations were carried out with the Wil- 
son action and using a relatively large value for the lattice 
spacing, a - 0.21 fm for SU(2) and a - 0.165 for SU(3). 
It is not clear yet whether this have any influence on such 
propagator. 

Recently, in [29j rigorous upper and lower bounds for 
the zero-momentum gluon propagator were derived and 
its scaling behaviour with the volume analyzed for the 
SU(2) Yang-Mills theory. The authors found a finite 
and non-vanishing D(0) for SU(2) in the infinite volume 
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limit. In [30(, the same analysis was performed but for 
the SU(3) Yang-Mills theory. It turns out that for SU(3), 
the infinite volume limit gives a D(0) = 0. The different 
behaviours in the infinite volume limit is a puzzling re- 
sult and, presently, we have no understanding for such a 
difference. 

In what concerns the lattice ghost propagator 0, [2l|, 
[22l . [25l . I31L [H, [3^, the simulations show an enhanced 
propagator but not in agreement with that predicted by 
the solution (T5|). The lattice data seems to be closer 
to the solution of [7|. Indeed, the lattice data shows a 
propagator which is almost identical to the perturbative 
propagator. 

In order to try to understand the difference between 
the DSE solution which compatible with the GZ and KO 
confinement mechanisms, in (40| the DSE were solved on 
a 4D symmetric torus. According to the authors, the 
gluon and ghost propagators approach slowly the infinite 
volume value. Moreover, they claim that to observe the 
suppression of the gluon propagator one should go to 
volumes as large a s (10 fm) 4 . Note that recent lattice 
simulations [2G|,[2a] have volumes well above the (10 fm) 4 . 

In (2lL I25H32I [38l [3^| we have tried to measure both the 
gluon and ghost propagators using a set of large asym- 
metric lattices, i.e. L 3 xT with T larger than L, to access 
the deep infrared region and to check the compatibility 
between the lattice data and the solution (J3j- ((4j) . The 
lattices used are larger than 10 fm, by a factor of ~ 2.5, 
in the temporal direction and are much shorter, by a fac- 
tor of ~ 1/5, in the spatial directions. For the L 3 x T 
lattices, besides the finite volume effects also observed in 
simulations with symmetric L 4 lattices, one has to care 
about how the asymmetry changes the gluon and ghost 
propagators data. The finite lattice effects in D(q 2 ) and 
G(q 2 ) in the asymmetric lattices are, qualitatively, equal 
to the effects observed in the solutions of the DSE on 
a symmetric 4D torus [4(| and on the simulation of 3D 
asymmetric SU(2) lattices [4l| . 

For the gluon propagator, in 39], using the lattice 
dressing function q 2 D(q 2 ) 1 and excluding the zero mo- 
mentum point from the analysis, we have demonstrated 
that the gluon lattice data is compatible with ([3]). We 
would like to call the reader attention that we started by 
simulating 16 3 x 128 and 16 3 x 256 lattices for (3 = 6.0. 
For these two lattices, it turned out that the 16 3 x 128 
data cannot be described by a pure power law. Only the 
fits to the 16 3 x 256 gluon data have acceptable \ 2 1 d.o.f '., 
i.e. x 2 /d.o./. < 2. Furthermore, an attempt to extrapo- 
late the lattice to infinite spatial volume suggests a K in 
the range 0.498 to 0.525. Despite this result, which gives 
some support to the KO and GZ confinement scenarios, 
before extrapolations the lattice data shows no suppres- 
sion of the gluon propagator for small momenta, except 
for D(0) when compared to the first non-zero momentum 
(see figure [T]). 

In what concerns the gho st propagator computed using 
asymmetric lattices pll |25|] . we have observed that G(q 2 ) 
is enhanced, compared to the perturbative solution, in 



TABLE I: Lattice setup. All simulations use a Monte Carlo 
sweep of 7 overrelaxation updates with 4 heat bath updates. 
The number of thermalization (Therm) and separation (Sep) 
sweeps refers to the combined sweeps. See [391] for details. 


Lattice 


Therm. 


Sep. 


# Conf. 


8 3 x 256 


1500 


1000 


80 


10 3 x 256 


1500 


1000 


80 


12 3 x 256 


1500 


1000 


80 


14 3 x 256 


3000 


1000 


128 


16 3 x 256 


3000 


1500 


155 


18 3 x 256 


2000 


1000 


150 



the infrared region (see figure [lj , but not as much as 
predicted by the DSE solution ©. 

In this work we discuss a method which, assuming a 
pure power law behavior for the infrared propagators, 
aims to measure the infrared exponents without relying 
on data extrapolation. The measure relies on the defini- 
tion of a convenient ratio of propagators which, in prin- 
ciple, suppresses the finite volume effects. The results re- 
ported show that, for the same asymmetric lattices used 
previously, the method provides estimates of the gluon 
propagator exponent which are stable against variation 
of the range of momenta and variation of the spatial lat- 
tice extent L. Furthermore, the measured k are compat- 
ible with the GZ confinement scenario, i.e. they predict 
a vanishng zero momentum gluon propagator. On the 
other hand, the data for the ghost propagator, although 
stable against variation of the lattice volume, is not com- 
patible with a pure power law behaviour. 

The paper is organized as follows. In section [ill we re- 
port on the asymmetric lattices used to investigate the 
infrared gluon and ghost propagators. In section [TTT1 the 
method used to measure the infrared exponents is dis- 
cussed together with the results for the gluon propagator. 
The section also includes a discussion of the finite volume 
and asymmetry effects. In section ITVl the method is ap- 
plied to the ghost propagator and finally in section fVl we 
draw the conclusions. 



II. DETAILS OF THE SIMULATIONS 

In this article we use Wilson action, /? = 6.0, gauge 
configurations for the lattices reported in table |TJ The 
asymmetric lattices have a temporal extension of T = 256 
(~ 26 fm) which allow us to access momenta as low as 
48 MeV. The main difference to [3!| being the larger 
statistics for the largest lattices. 

The propagators were computed in the minimal Lan- 
dau gauge and the gauge fixing was performed using a 
Fourier accelerated steepest descent algorithm; see [39| 
for details and definitions. 

In the following, we will discuss the infrared propaga- 
tors and will consider only time-like momenta, defined 
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FIG. 1: Bare gluon and ghost propagators for time like mo- 
menta. Note the logarithmic scale for the ghost propagator. 
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2 / nn\ 
q[n\ = q 4 [nj = - sin ^— j , n = 0, 1, 

where T is the time lattice extent. For the conversion to 
physical units we use a" 1 = 1.943(47) GeV [12]. 

The ghost pro pagator was computed with the method 
described in [43|], for the smallest q[n]. In the calcula- 
tion of the D(q 2 ) and G(q 2 ), the statistical errors were 
evaluated with the jackknife method. Otherwise, the sta- 
tistical errors were computed using the bootstrap method 
with a 68% confidence level. The bare lattice gluon and 
ghost propagators are reported in figure [TJ 



III. THE GLUON PROPAGATOR 

For the measurement of the infrared exponents, it will 
be assumed that the lattice dressing functions Z(q 2 ) and 
F(q 2 ) are described by pure power laws {q 2 ) a , as in © 
and ((4]), times a factor A(g) which summarises the fi- 
nite volume corrections and/or deviations from the pure 
power law. If these corrections are constant (small), they 
are eliminated (suppressed) by taking ratios of Z and F 
at consecutive lattice momenta, i.e. 



hi 



G(q 2 



1]) 



Q{q 2 [n]) 



= a In 



q 2 [n] 



(6) 



and, according to ([3]) , 
2^ 



where for the gluon Q{q 2 ) = Z(q 2 

a = 2k. For the ghost G(q 2 ) = F(q 2 ) and a = — k 
In the infrared one expects that corrections to © are 
subleading. Defining 



R z [n] = In 
R q [n] = In 



Z(q 2 [n + 1}) 



Z{q 2 [n}) 
q 2 [n + l] 
q 2 [n] 



(7) 
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FIG. 2: Ratios of gluon (left) and ghost (right) dressing func- 
tions. The statistical errors were computed using 1000 boot- 
strap samples for 8 3 — 12 3 x 256 lattices and 1800 bootstrap 



samples for 14" - 18 J x 256 lattices. 

we get, for the gluon propagator, 
R z [n] = 2k R q 



(8) 



The gluon data for Rz [n] as a function of R q [n] , see fig- 
ure shows a linear behavior for a surprising large range 
of momenta for all the lattices 8 3 — 18 3 x 256. Moreover, 
the slopes seem to be similar for all lattices. It seems that 
the corrections to ([3]) not suppressed by the ratios show 
up as a constant which adds up to ©. This hypothesis 
can be tested fitting the ratios to 



R z [n] = 2KR q [n] 



C 



(9) 



assuming that C is a constant. The fits of the lattice ra- 
tios to © using momenta up to ~ 400 MeV are reported 
in table HH 

The measured k's are stable against variation of the 
fitting range and spatial lattice size. Note that the k 
measured with the smallest lattice 8 3 x 256 is compatible 
with the k measured using our largest volume 18 3 x 256. 
Moreover, the central values for k are clearly above 0.5 
and typically, within one standard deviation, k > 0.5. 
In particular, for the three largest volumes 14 3 — 18 3 x 
256, which have the largest statistics, k > 0.5 within one 
standard deviation. Statistical errors on k decrease as 
the fitting range increases, because one is using a larger 
set of data. For the largest fitting range (q < 381 MeV) 
and for the three largest lattices, the k are compatible 
with 0.5 only within 4 standard deviatons. In this sense, 
in what concerns the infrared gluon propagator, the fits 
to ([9]) point towards k ~ 0.53, suggesting a vanishing 
gluon propagator at zero momentum. 

Note that the absolute value of the constant C, in gen- 
eral, approaches zero as the lattice volume increases (see 




1/V R [n] 



FIG. 3: C as function of 1/V for asymmetric lattices with T = 
256. Values were taken from table [TT] considering the largest 
fitting range q < 381 MeV and the smallest q < 191 MeV. 



figure [3]) , opening the possibility of having a vanishing C 
in the infinite volume limit. In this sense, C can be seen 
as a measure of the finite volume corrections to the pure 
power law. 

In conclusion, the results reported in this section give 
support to an infrared gluon propagator behaving as a 
pure power law, with 

^ 2 )-('7 2 ) 2 ^ 1 -(<Z 2 ) 06 . (10) 

Accordingly, this means a vanishing zero momentum 
gluon propagator in agreement with the GZ confinement 
mechanism. 



A. Finite Volume and Asymmetry Effects 

So far we have used only data from asymmetric lattice 
simulations with T = 256. As a check to the above re- 
sults, in this subsection we investigate the finite volume 
effects and the changes introduced from using a time-like 
direction much larger than the spatial directions. There- 
fore, besides the asymmetric lattices of table HI we are 
going to consider the data coming from the lattices re- 
ported in table Hill 

For on-axis momenta, the gluon propagator data (not 
shown here) for the lattices reported in tables U and IIIII 
show very clearly finite volume and asymmetry effects. 
However, here we are interested in the ratios Rz as a 
function of R q , which are reported in figure [3J 

Again, the data shows a linear behaviour for all the 
lattices of table IIIII A fit of the ratios computed using 
the 16 3 x 128 data to @ for the smallest fitting range 
(q < 381 MeV) gives K = 0.541(19) and C = -0.239(38) 
with a x 2 /dof = 0.01. The k value is in excellent agree- 
ment with the corresponding figure for 16 3 x 256. For the 
C figures, it comes that C128 — 2 x 6256- This result re- 



FIG. 4: Rz as function of R q for the lattices of table IIIII 
and the two largest asymmetric lattices of table [I] The figure 
contains data for q < 500 MeV for 16 3 x 256, q < 900 MeV 
for 16 3 x 128 , q < 1.5 GeV for 48 4 and = 6.0, q < 1.2 GeV 
for 64 4 and f3 = 6.0 and q < 1.6 GeV for 64 4 and (3 = 6.2. 



inforces our interpretation of C as a measure of the finite 
volume corrections to the pure power law. 

In what concerns the symmetric lattice data in fig- 
ure 21 despite their relative large smallest momentum, 
in the infrared the slope follows essentially the slope of 
the asymmetric ratios. The reader should be aware that, 
due to the q m in of the symmetric lattices considered here, 
extracting infrared exponents by fitting the ratios for 
the symmetric lattices to equation (JSJ) could be some- 
what meaningless. For each of the symmetric lattices, 
one can estimate a k directly from the two most right 
points which gives k ~ 0.6 or higher. The main difference 
seems to be in the constant C, which takes larger values 
for smaller lattices. In particular, note that (3 = 6.2, 
64 4 and (3 — 6.0, 48 4 have essentially the same physical 
volume and their ratio data is, within statistical errors, 
indistinguishable. 

Given the results summarized in figure 31 one can con- 
clude that the ratios defined in ([7]) provide reliable esti- 
mates of the gluon infrared exponent, i.e. the estimates 
seem to be independent of the lattice volume and asym- 
metry. 



IV. THE GHOST PROPAGATOR 

In what concerns the ghost propagator, the data for the 
ratios of the dressing functions, see figure [2 do not show 
a linear behaviour as in the case of the gluon propagator. 
We have tried a number of functional forms to fit the 
data, but their \ 2 /d.o.f. was always too large. At most, 
the slope of ratios of the ghost dressing function suggests 
a negative value for the ghost infrared exponent. 

In figure [5] we show the ratios of the ghost dress- 
ing functions for the three larger lattices, including the 
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TABLE II: Fitting the gluon ratios with equation © for L 3 x 256 lattices. The first line is the maximum momentum used in 
the fit. \ 2 stands for \ 2 /d.o.f.. The errors in k are statistical and were computed with the bootstrap method. 





Qmax • 


1Q1 A/TeV 




zoo ivie v 




OCR MeV 

zoo ivie v 




000 ivie v 




ool ivie v 




L 




Par am. 


X 


Par am. 


X 


Par am. 


X 


Par am. 


X 


Par am. 


,,2 

X 


8 


K 


0.526(27) 


0.12 


0.531(19) 


0.11 


0.531(13) 


0.08 


0.522(16) 


0.48 


0.527(12) 


0.54 




C 


-0.179(54) 




-0.194(34) 




-0.193(19) 




-0.171(28) 




-0.184(18) 




10 


hi 


0.511(35) 


0.69 


0.531(25) 


0.98 


0.525(21) 


0.74 


0.523(17) 


0.56 


0.527(16) 


0.50 




c 


-0.114(66) 




-0.161(42) 




-0.146(30) 




-0.144(21) 




-0.150(19) 




12 


hi 


0.509(31) 


0.11 


0.517(21) 


0.16 


0.508(18) 


0.33 


0.521(18) 


0.84 


0.530(14) 


1.03 




c 


-0.094(56) 




-0.112(35) 




—0.094(25) 




-0.119(27) 




—0.138(18) 




14 


hi 


0.536(24) 


0.33 


0.540(19) 


0.20 


0.548(16) 


0.39 


0.545(12) 


0.34 


0.542(11) 


0.34 




c 


-0.114(44) 




-0.123(30) 




-0.140(21) 




-0.134(15) 




-0.127(12) 




16 


K 


0.539(22) 


1.77 


0.528(17) 


1.24 


0.534(12) 


0.96 


0.536(12) 


0.78 


0.539(11) 


0.68 




c 


-0.125(43) 




-0.102(30) 




-0.112(19) 




-0.118(14) 




-0.123(12) 




18 


K 


0.529(20) 


0.39 


0.516(16) 


0.77 


0.523(14) 


0.85 


0.536(11) 


1.79 


0.5398(95) 


1.58 




c 


-0.099(36) 




-0.068(25) 




-0.085(19) 




-0.111(14) 




-0.119(13) 





TABLE III: Lattice setup used to check finite volume and 
asymmetry effects. To convert to physical units we used the 
central values for the lattice spacing reported in [42| : a = 
0.1016 fm for (3 = 6.0 and a = 0.0728 fm for f3 = 6.2. Volume 
is reported in fm. q m in (in MeV) is the smallest non- vanishing 
momentum for each lattice. For symmetric lattices only on- 
axis momenta are considered. 





Lat 


Volume 


Qmin 


# Conf. 


6.0 


16 3 x 128 


(1.63) 3 x (13.00) 


95 


164 


6.0 


48 4 


(4.88) 4 


254 


104 


6.0 


64 4 


(6.50) 4 


190 


120 


6.2 


64 4 


(4.66) 4 


266 


99 



R F [n] = - nR q [n] + C 



(11) 



where k = 0.529 and C was adjusted to reproduce the 
ratio computed using our smallest momenta. The figure 
shows that either the data is still far from the linear be- 
haviour or the infrared ghost propagator does not follow 
a pure power law. 



V. CONCLUSIONS 

In conclusion, in this article we discuss a method which 
does not rely on extrapolations to the infinite volume and 
i) allows to check if the infrared gluon and ghost propa- 
gators follow a pure power law; ii) measure the infrared 
exponents if the propagator is compatible with a power 
law. 

For the gluon propagator the lattice data validates the 
power law behavior and the corresponding infrared ex- 
ponent is measured. From the two large asymmetry lat- 
tices and the two smallest fitting ranges we have k = 
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FIG. 5: -RfM as a function of -R q [n] for the three largest 
lattices. The dash lines show, for each lattice, the curve 
— hiR q [n] + C where k = 0.529 and C adjusted to reproduce 
the right end point in the graph. 



0.529(20) and 0.539(22) which gives a combined value of 
k = 0.535(14). As discussed, the ratio method devised in 
this paper seems to be independent of the asymmetry and 
lattice volume, at least for the simulations reported here. 
This gives confidence on our estimate for k. Furthermore, 
the reported exponent k = 0.535(14) is within the range 
of values allowed by the renormalization group analysis 
of OH!]. Given the results for the ghost propagator, it 
is not clear to the authors if this k should be or not in 
that range. More, the measured k supports a vanishing 
zero momentum gluon propagator in agreement with the 
Gribov-Zwanzigcr gluon confinement scenario. 

The reader should note that, on a finite lattice, simula- 
tions have produced always a finite and non-zero value for 
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the gluon propagator at zero momentum. The method 
discussed here do not use directly -D(O). Anyway, as men- 
tioned before, the conclusions are in perfect agreement 
with the scaling analysis of -D(O) as a function of the 
volume performed in [30j], which gives a D(0) = in the 
infinite volume limit. In this sense, the results of [3(3] rec- 
onciliate the results of the lattice simulations with those 
reported in this work, i.e. that in the Landau gauge the 
continuum result is D(0) — 0. 

In what concerns the ghost propagator, the lattice data 
suggests that the propagator does not follow a pure power 
law. In this sense, our simulation does not validate the 
infrared solution ([4]) of the DSE. The observed difference 
between the lattice and the DSE solution could be be- 
cause either the lattices used here are not long enough or 
the ghost propagator does not follow a pure power law 
in the deep infrared region. The answer to this question 
requires further investigations. 
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APPENDIX A: FINITE VOLUME EFFECTS ON 
THE INFRARED GLUON DRESSING 
FUNCTION: AN EFFECTIVE 
PARAMETERIZATION 

The results on the ratios of the gluon dressing function, 
discussed in section IIII1 suggest that the finite volume 
and asymmetry corrections are summarized in C . Let 
A{q) be the multiplicative correction to the continuum 
dressing function Z(q 2 ), i.e let us assume that the lattice 
dressing function is given by 

Z Lat (q 2 ) = Z(q 2 )A(q). (Al) 

Then, A(q[n + 1]) = A(q[n])e which allows to write 

dA(q) ^ A(q[n + 1])-A(q[n}) 
dq q[n + 1] - q[n] 

~ A(q)^- = A(q) A (A2) 
aT 

where A is a constant. The integration of the last equa- 
tion gives A(q) = Aq exp(Aq) where Aq is a constant 



of integration. The reader should note that this form for 
A(q) is only valid in the infrared region. Computing A(p) 
for the full range of momenta is beyond the scope of this 
work. Nevertheless, for high momenta, where the propa- 
gator is well described by its perturbative behavior and 
the lattice gluon propagator does not seem to depend on 
the volume, one should have A(p) ~ 1. In principle, this 
restriction allow us to fix Ao via the normalization condi- 
tion A(gh) = 1 where qh is a sufficiently high momentum 
value. The present lack of knowledge on the functional 
form for A(p) at momenta q > 400 MeV prevent us to 
compute Ao. At best, one can estimate Ao at the highest 
momentum q (~ 400 MeV) where one can still approx- 
imate A(q) = Aq exp(Aq). Then, A = exp(— Aqh) 
and 

A(q) = A exp(Aq) ~ exp (A(q - q h )) (A3) 

for the infrared region. 

For practical purposes, like fitting the lattice data, A 
can be absorbed into the definition of u>, and we get an 
exponential correction to Z(q 2 ), 

Z Lat (q 2 ) = Uj(q 2 ) 2K e A \ (A4) 

with the constant A parametrizing the finite volume ef- 
fects and/or the corrections to the pure power law and k 
is the continuum exponent. The results of fitting (|A4|) to 
the lattice gluon dressing function are reported in table 

El 

The k values in tables |TT] and IIVI are essentially the 
same. The good agreement between the two results can 
be viewed as a cross-check of the original method dis- 
cussed previously and, in this sense, gives further confi- 
dence in the final result. The reader should note that A, 
like C, takes negative values and approaches zero as the 
lattice volume increases. 

In what concerns the function A(q), see eq. (|A3|) . 
it provides a correction to the gluon propagator which 
decreases with increasing momenta, as expected if it 
models the finite volume effects. Indeed, A(q = 0) ~ 
cxp(— Aqh) » 1 and A{qh) ~ 1 - see figure M where 
we have corrected the propagator using qh — 400MeV. 
The reader should not forget that we are only estimating 
the corrections due to finite volume /asymmetry. Figure[5] 
shows a corrected gluon propagator which is supressed in 
the infrared region. Moreover, the corrected data shows 
also a good agreement between all propagators, excluding 
our smallest lattice, for momenta below 200 MeV. 
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TABLE IV: Fitting the gluon dressing functions to (|A4|I for L 3 x 256 lattices. The first line is the maximum momentum used 
in the fit. \ 2 stands for \ 2 /d.o.f.. The errors are statistical and were computed with the bootstrap method. A in GeV -1 . 
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FIG. 6: Corrected and uncorrected gluon propagator using 
eq. (0 with q h = 400 MeV. 
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